Universal scaling exponents in shell models of turbulence: 
Viscous effects are finite-size corrections to scaling 
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In a series of recent works it was proposed that shell models of turbulence exhibit inertial range 
scaling exponents that depend on the nature of the dissipative mechanism. If true, and if one could 
imply a similar phenomenon to Navier-Stokes turbulence, this finding would cast strong doubts on 
the universality of scaling in turbulence. In this Letter we propose that these "nonuniversalities" 
are just corrections to scaling that disappear when the Reynolds number goes to infinity. 



PACS numbers 47.27.Gs, 47.27.Jv, 05.40.+j 

The aim of this Letter is to question a growing con- 
sensus concerning the non-universality of the scal- 
ing exponents that characterize correlation functions that 
appear in the context of numerical studies of shell mod- 
els of turbulence. Shell models, and in particular the 
so-called GOY model H|), became very popular due to 
their ease of simulation and the fact that they appear 
to exhibit anomalous scaling that is strongly reminiscent 
of the behavior of the scaling exponents of Navier-Stokes 
turbulence @-||. The GOY model is a simplified reduced 
wavenumber analog to the spectral Navier-Stokes equa- 
tions. The wavenumbers are represented as shells, each of 
which is defined by a (real) wavenumber k n = fcoA™ where 
A is the "shell spacing" . There are N degrees of freedom 
where N is the number of shells. The model specifies 
the dynamics of the "velocity" u n which is considered a 



complex number. 

- d 
-dt 

k n 
~ ~2 



, N. The GOY model reads 



= fn + i[k n +lU n +lU n +2 



-1U„ 



-U n -lU n -2\ 



(1) 



In this equation * stands for complex conjugation, v is 
the "viscosity" , and to motivate the discussion we already 
specified a one parameter family of models in which the 
exponent a determines the type of viscous dissipation. 
In the following we use A = 2 and a forcing restricted to 
the first two shells. 

The scaling exponents can be defined using the "corre- 
lation functions" (|u„| p ). It was shown however that the 
GOY model suffers from periodic oscillations superposed 
on the scaling laws obeyed by these functions jl(J . It was 
proposed in Q that it is advantageous therefore to focus 
on the scaling properties of correlation functions of the 
type 



F q (k n ) 



Im I u n u n+1 u n+2 + -w„_iw„u„ + i 



These functions scale with fc„, F q (k n ) ~ kn^ as long as F q as a function of k n , for different values of the param- 
k n is smaller than some dissipative cutoff scale kd and suf- eter a = 1, 2, 3 and N — 22. The dissipative cutoff 



9/3 



(2) 



ficiently larger than k\ . The exponents C, q are the scaling 
exponents whose universality is questioned. The defini- 
tion (||) has the advantage that £3 is known exactly (if 
one assumes that inertial range scaling is universal) and 
it is Cs = 1 I- 

The issue to be discussed in this Letter is introduced in 
Fig. 1 which shows the results of numerical simulations 
of this family of models. The simulations were performed 
using the adaptive step-size backward differentiation al- 
gorithm (DDEBDF) from the SLATEC library 011 ' 
We imposed a random forcing on the first two shells, 
fi= 5 10~ 3 and f% = af\ where a typical value of a is 
0.75. The time correlation of the forcing has been chosen 
to be of the magnitude of the natural forcing turn-over 
time-scale r = \j\Jk\f = 40. 
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FIG. 1. The correlation functions Fi to F5 averaged over 
about 1500 forcing turn-over scales for a model of 22 shells 
with a dissipative term proportional to k 2 , (y = 5 X 10 -7 ), 
fc 4 , (y = KT 14 ) and fc 6 , (y = 2 x 10" 22 ). 



We plot in log-log coordinates the correlation functions 
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Table 1. Simple scaling regression k n for regular and 
hyperviscosity, a — 1,2,3. 
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22 
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0.69 


0.97 
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fc 6 


22 


0.34 


0.62 


0.85 


1.06 


1.25 



is chosen (by fixing the parameter v) to keep kd « ko2 17 
for all values of a. In Table 1 we show the results ob- 
tained by simple linear fits between the shells 4 and 13. 
Results of this type were interpreted in previous works 
as an indication that the inertial range scaling 
depends on a, casting therefore severe doubts on the 
universality of the scaling exponents Cq- It should be 
stressed that this is not a trivial issue; if indeed the scal- 
ing exponents depend on the dissipative mechanism this 
could cither imply a major departure from the standard 
thinking about universality in turbulence, or an indica- 
tion that shell models are fundamentally different from 
Navier-Stokes dynamics. The main point of this Letter is 
that none of this is necessary. We propose that these scal- 
ing plots should be interpreted as standard corrections to 
scaling that can be factored away by taking into acount 
the effect of the viscous dissipation. These effect becomes 
less important when the number of shells increases and 
when the length of the inertial range increases, as we 
show below in detail. 

To provide strong evidence that we are faced with finite 
size effects we integrated the equations of motion for dif- 
ferent number of shells. We chose three sets of data with 
TV = 22, 28 and 34 shells. In all cases the integration was 
performed over 1500 forcing turn-over time scales. The 
viscosities were chosen so that the width of the dissipative 
range is about 5-6 shells: for a = 1 (normal dissipation), 
v = 5 x 10" 7 , 8 x 10~ 9 and 4 x 10~ n for N = 22, 28 
and 34 respectively; for a = 2, v = 10~ 14 , 2 x 10~ 20 and 
2 x 10~ 26 and for a = 3, v = 2 x 10~ 22 , 10~ 30 and 10~ 41 . 
The results of this simulation are a direct illustration of 
the finite size effect: for any value of a, we can replot 
our data for N = 22, 28 and 34, moving the abscissa to 
coalesce the near dissipative region. Consider for exam- 
ple Fs(k n ) for which we have a theoretical expectation, 
i.e., F^,{k n ) ~ k~ l in the inertial range, but for which 
the data in Table 1 indicate significant deviation from 
£3 = 1. In Fig. 2 we present double-logarithmic plots of 
k n F^{k n ) vs. k n for a = 2. We expect this function to be 
constant in the inertial range. Instead, we see a strong 
viscous effect, leading to deviation from constancy over 
10 shells even before deep dissipation sets in. Increas- 
ing the number of shells moves the problematic region to 
higher shells. We have shifted the abscissa in Fig. 2 to 
collapse all our data together. The data indeed collapses, 
with an N independent (fixed size) cross over region. 
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FIG. 2. Log-log plots of k n Szikn) vs. k n in case of hyper- 
viscosity of index a = 2 with different numbers of shells and 
viscosities. The collapse has been obtained by shifting the ab- 
scissa. In the inner frame, we zoom on the y-axis to exhibit 
the finite size effect occuring in the inertial range. The solid 
line shows the constant behavior expected theroretically. One 
observes clearly that the departure from this constant value 
only occurs in a region of about ten shells near to the vis- 
cous transition. When the inertial range is large enough, the 
predicted behavior is recovered. 

This is a sure indication that we are dealing with a finite 
size effect and not the breaking of universality. Any devi- 
ation of the apparent (3 from unity is going to disappear 
slowly when the size of the inertial range increases. Of 
course, the linear fit procedure is very sensitive to the 
existence of a bump of the type shown in Fig. 2 and 
the determination of the exponents in Table 1 is very 
unreliable. 

As is commonly found in scaling systems that suffer 
strong finite-size effects, it is advantageous to fit the whole 
scaling function rather than its power law part alone. In 
other words, to extract reliable inertial range scaling ex- 
ponents, we propose to fit the correlation functions over 
their full domain, including the deep viscous regime. To 
achieve such a fit, we need to consider first the viscous 
range, where we may guess a generalized exponential de- 
cay: 



exp 



(3) 



where kd is an appropriate dissipative scale. In the sta- 
tionary state we expect a balance between the non-linar 
transfer and the dissipatve damping for every shell. Due 
to the exponential decay, the most relevant non-linear 
term for the n'th shell is fc„_iM„_iU n _2, which has to 
balance with fc 2a u„. Up to logarithmic corrections, the 
exponent x is then determined from 
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FIG. 3. Plots of log 2 | log 2 Ss(k n )\ vs log 2 (fc„) for indices of 



hyperviscosity a = 1,2,3 and viscosity u = 5 x 10 



10" 



and 8 x 10~ . The dotted lines represent the asymptotic 
stretched exponential behavior. The slopes are respectively 
0.69, 0.76 and 0.86 to compare to the theoretical expectation 
x = 0.694. 
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FIG. 5. Third order exponent £3 computed via the fitting 
procedure described in the text for a = 2 and 3 obtained 
for N = 22, 28, 34. The results have been plotted against 
l/(An) 2 where An is the estimated width of the inertial 
range. The straight lines (with respective slope -3.5 and -14) 
indicate that the calculated exponents converge linearly to- 
ward the predicted value £3 = 1. 



PI 



60 

o 



0.0 



-20.0 



-60.0 



-100.0 



-120.0 



-140.0 



oF 3 (o=l) 

Fit-Fj (a=l) 

□ F, (a=2) 

1 it 1 .« :•■ 

o F 3 (a=3) 
Fit-F, (a=3) 




-5.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0 

l°g 2 (k) 

FIG. 4. Third order structure function Fz(k n ) for a = 1, 2 
and 3. The symbols indicate the results of simulations with 
34 shells; the lines are the fits performed between the 3rd and 
32nd shells. The results have been shifted in the y-coordinate 
for clarity. 
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or equivalcntly, 



1 + \ x - A 



x \2x 



so that we have finally 

x = log A 



1 + V5 



(4) 
(5) 

(6) 
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FIG. 6. Same figure as above for the second order expo- 
nent £2- The illustrative lines have slopes -6 and -12. The 
final value ("2 = 0.711 has been obtained in the normal dissi- 
pation case a = 1. 



Our simulations provide excellent support for the deep 
dissipative stretched exponential dependence, but the 
predicted value of x depends on the coefficient a, see 
Fig. 3. In the inertial range we expect a power law, and 
to estimate the behavior in the cross over region we note 
that the ratio of the viscous to nonlinear term can be 
written in the form 



3 



vk 



2n „ 



vk 



2a. 



k d u 2 d 



vk 2 fud k r , 



2Q-2/3 



(7) 



where Ud is the velocity at the viscous cross-over shell 
with k n — kd- For simplicity we assumed u n ~ k n ^ 3 
for this estimate. Taking (Q) as a first order correction 
of the inertial range power law by viscous effects we can 
propose a functional form for F q (k n ) over its full domain: 



F q (k n 



2a-2/3 



exp 



(8) 



By construction this form fits the scaling law and the 
deep dissipative form for k n <<C kd and k n ^> kd re- 
spectively. For crossover values of k n the formula offers 
two additional free parameters, i.e. f3 q , kd, q - This is the 
minimal number of fit parameters that is needed to in- 
terpolate between the inertial and the dissipative ranges. 

As an example of the excellent fit that this formula 
provides we show in Fig. 4 F^(k n ) for a = 1, 2 and 3. 
This simulation employed 34 shells, and the fit were per- 
formed between the 3rd and 32nd shells. Similarly good 
fits are found for all the data that we analyzed. 

Accordingly, we recorded the apparent values of the 
scaling exponent £ g for increasing values of the N, and 
noticed that they converge to a given value. This conver- 
gence can be seen in Figs. 5 and 6 in which the scaling 
exponents (2 and £3 are plotted for a = 2, 3 as a func- 
tion of l/(An) 2 where An = log 2 (fcd/fe) is the estimated 
width of the inertial range, and kd is obtained from the 
fit. It is obvious that the exponents converge towards a 
universal value that agrees with a = 1. 

The conclusion of this Letter is that it is dangerous 
to measure scaling exponents in multiscaling situations 
without taking into account crossover and finite size ef- 
fects. Even in the case of shell models which appear to 
have much longer inertial ranges than Navier-Stokes tur- 
bulence, simple log-log plots arc misleading. We propose 



that in future determination of scaling exponents in nu- 
merical experiments similar care must be given to such 
issues. 
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